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Abstract 



Janus and Epimetheus are famously known for their dis- 
tinctive horseshoe-shaped orbits resulting from a 1:1 or- 
bital resonance. Every four years these two satellites swap 
their orbits by a few tens of ki lometers as a re s ult o f 
their close encounter. Recently Tiscareno et al. ( 2009( ) 



have proposed a model of rotation based on images from 
the Cassini orbiter. These authors inferred the amplitude 
of rotational librational motion in longitude at the orbital 
period by fitting a shape model to the recent Cassini ISS 
images. By a quasiperiodic approximation of the orbital 
motion, we describe how the orbital swap impacts the 
rotation of the satellites. To that purpose, we have de- 
veloped a formalism based on quasi-periodic series with 
long and short-period librations. In this framework, the 
amplitude of the libration at the orbital period is found 
proportional to a term accounting for the orbital swap. 
We checked the analytical quasi-periodic development by 
performing a numerical simulation and find both results 
in good agreement. To complete this study, the results 
regarding the short-period librations are studied with the 
help of an adiabatic-like approach. 

1 Introduction 

The orbital motion of Janus and Epimetheus presents 
a peculiar horseshoe-shaped orb it resulting from a 

1:1 o rbital resoiiance (e.g; Dermott and Murray 

19811: lYoder et all Il983l ; iMurrav and DermottI llQQol : 
Jacobson et al.l 20081 and references therein) . Every four 
years the two satellites swap their orbits by a few tens of 
kilometers as a result of their close encounter. As the mass 
of Janus is 3.6 times greater than the mass of Epimetheus, 
the dynamical motion of the latter is more sensitive to the 



swap than the dynamical motion of Janus. 

The rotational motion of the satellites depends mainly 
on the gravitational torque of Saturn acting on the dy- 
namical figure of each moon. The expression of the grav- 
itational torque is: 



1 ~ ^ — u X [l\u 



(1) 



with G the gravitational constant, Ms the mass of Saturn, 
[/] the inertia tensor of the moon, r the distance between 
Saturn and the moon, and u the unit vector toward Saturn 
in the moon's reference frame. The gravitational torque T 
depends on the relative Saturn-moon distance, hence the 
swap also yields his signature on the rotational motion of 
the satellites. 

First estimates of the rotational motion of the two coor- 
bit al satellites Janu s and E pimetheus have been obtained 
by iTiscareno eV^ (|2009[) . From images provided by 
the Cassini orbiter, they fitted a numerical shape model 
of the moons, which included the amplitude of the li- 
bration in longitude. The libration in longitude corre- 
sponds to_the_oscinationof_^ body along its equatorial 
plane. Tiscareno et al. ( 20091 ) obtained an amplitude of 
5.9° ± 1.2° for Epimetheus. For Janus, the uncertainty 
on the fit of the libration determination is too large to 
yield an accurate librational amplitude. However, based 
on their shape model, Tiscareno et al. suggested a value 
of 0.33° ±0.06° for the amplitude of the libration in longi- 
tude. In addition, they identified an unexplained constant 
phase of 5.3° ± 1° for Janus, whereas for Epimetheus such 
offset is in the error bar. A recent numerical study by 
Novelles ( 2010[ ) explored the three-dimensional rotational 
motion o f these satellites based o n the numerical shape de- 
duced by|Tiscareno^£^ (2009). Noyelles' study suggests 
a strong influence of the swap on the rotational motion 
of Janus and Epimetheus, which we propose to explore in 
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the present study. 

The orbital motion of Janus and Epimetheus appears 
to be very regular, at least over a timescale of several 
thousands of years. Thus, the trajectories of these satel- 
lites can be considered as quasi-periodic. Schematically 
these trajectories evolve on three different timescales. The 
shortest corresponds to the mean motion of the satellites 
with a period of about 0.7 days. The second component 
has a period of 8 years and is associated with the close 
encounters of the satellites. The long-period component 
is the secular variations of the satellites eccentricities and 
inclinations over periods of a few thousands years. Since 
the rotation of Janus and Epimetheus is synchronous with 
their orbital motion, it reflects these different timescales. 

One of the goals of this paper is to understand the in- 
fluence of the 8-year horseshoe motion on the rotational 
librations of Janus and Epimetheus. To this purpose, 
we develop an analytical solution of the rotation of these 
bodies that can simulate the main features of their spin. 
Then, in Section |3l we model the orbits of the two satel- 
lites through quasi-periodic expansions. The fourth sec- 
tion is dedicated to the description of the librational mo- 
tion of each satellite. We develop three approaches to 
describe the librational motion in details: (1) a quasi- 
periodic development that highlight the fundamental fre- 
quencies involved in this problem, (2) an adiabatic in- 
variant approach that focuses on the short-period libra- 
tions, and (3) a numerical approach to reach high accu- 
racy. Then we discuss the offset of Janus's orientation 
and investigate the effect of high spherical harmonics (or- 
der 3) and tidal coupling. We also discuss the influence 
of triaxiality on libration amplitudes. 



2 Physical librations for Keple- 
rian orbit 

First, let us recall the librational response of a satellite in 
synchronous spin-orbit resonance with a fixed Keplcrian 
orbit (constant semi- major axis a and eccentricity e). The 
position of the moon is determined through its relative 
distance with Saturn r, and its orbital longitude is defined 
by the draconic true longitude (angle between the body 
and its line of nodes) denoted by v. The mean anomaly 
will be denoted by £, while the true anomaly is denoted by 
/. We neglect the effect of the obliquity, which is small 
( Novellej . l20inh . so that the orientation of the body is 
specified by the angle 9 defined with respect to the line of 
nodes. 

The dynamical equation governing the rotation of the 
moon is the angular momentum balance with the grav- 
itational torque exerted by Saturn. If we introduce the 



physical libration 7 by the relation 9 ~ £ + uj + where 
the angle 7 represents the oscillations around the uniform 
synchronized rotation and w the fixed argument of the 
pericenter, the angular momentum equation reads: 



7+y(^) sin2(7 + ^-/) = 0, 



(2) 



The frequency a is the frequency of the free libration, 
also called frequency of the proper libation. It is equal 
to cr = n^3{B - A)/C where A < B < C are the nor- 
malized moments of inertia of the satellite and n its mean 
motion. 

For a small eccentricity e, the difference between the 
true and the mean anomalies is approximated at first- 
order in e by / — ^ = 2e sin £. In addition, we approximate 
a = r, and, for 7 small, the linearized equation of Eq. ([2]) 
is 

7 4-cr2^ = 2ecr2sin£ ^3-) 
The librational solution is then simply 



7 = sin [at + 1 



2ea2 



(4) 



where A^ , (j)-^ depend on the initial conditions and the 
right-hand side term is the forced libration. The forced 
libration oscillates at the mean motion frequency, and its 
amplitude is proportional to the ratio of the eccentric- 
ity to the difference between the square of the libration 
frequency a and the square of the forced frequency n. 
Therefore, the amplitude of the forced libration depends 
on both the magnitude of the torque and the proximity 
of the free libration frequency to the orbital frequency. 

In the following sections, we describe the orbits and 
investigate the impact of the horseshoe-shaped orbit on 
the physical librations. 

3 Orbital description of Janus and 
Epimetheus 

3.1 Osculating elliptical elements and 
fundamental frequencies 

The co-orbital satellites Janus and Epimetheus are fa- 
mously known to exchange their orbits every four years. 
This swap takes a short time-span, which does not exceed 
six months. In order to model the peculiar orbital motion 
of these satellites, we numerically integrate the three-body 
problem composed of Saturn, Janus, and Epimetheus, in- 
cluding the oblateness J2 of Saturn By using the fre- 
quency analysis developed by iLaskaii (|l988l 120051 ) for the 
purpose of Celestial Mechanics studies, we express the nu- 
merical solution as a quasi-periodic function expanded in 
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Fourier scries, where each frequency is a hnear combina- 
tion (with integer coefficients) of six fundamental frequen- 
cies (proper frequencies) denoted by {n,v,gj,gE,sj,SE) 
(in the following, we use the subscript J for Janus and E 
for Epimetheus). The first of these frequencies, n, called 
proper mean motion, is associated to the mean orbital 
motion common to the two co-orbital satellites, and it is 
constant along satellite orbits. In the same way, the five 
other proper frequencies are also con stant a,nd mig ht be 
considered as integrals of movement ( Laskarl . 120051 ) . The 



second fundamental frequency v corresponds to the libra- 
tion frequency along the horseshoe orbits, while the four 
last ones are associated to the motion of the pericenters 
{gE,gj) and of the ascending nodes {se,sj). 

If the gravitational interactions between Janus and 
Epimetheus were negligible, then the precessions of the 
pericenters and of the nodes, which arc only due to J2, 
would satisfy the relations gj = Qe and sj = se, and the 
eccentricities and inclinations of the two satellites would 
be constant (excepted for small short-period variations). 
However, due to the satellites mutual interactions, the 
proper frequencies of precessions of the orbits are slightly 
different, and both eccentricities and inclinations undergo 
large long-term variations whose periods are about 2760 
and 7690 years, respectively. The long-term variations of 
the eccentricities will be neglected in the present work. 
Their impact on the rotation of the co-orbital satellites is 
mentioned in section W?^ 

The initial conditions of the numerical integration s 



come from the ephemeris Horizons (jGiorgini et al 
The model used in the ephemeris Horizons takes into 
account the gravitational interactions with the Sun 
and the other main s atellites of the Saturnian system 
2008h . Thus, in order to be consistent 



( Jacobson et al. 



with that model, we fit the initial conditions of our inte- 
gration in such a way that the fundamental frequencies 
n and v arc the same in both cases. The entire set of 
fundamental frequencies, as well as their associated pe- 
riods, are displayed in Table [T] The frequencies ((?j,5_e) 
and {sj^se) are very closed indicating that the mutual 
interactions are small. 

Due to the high orbital precession rate generated by 
Saturn's oblateness, we define the elliptic elements by in- 
troducing Saturn's J2 in the third Kepler's law, which 
yields: 



n a = /i, with 

^i = G{Ms + mj + rns) il + ^jj^ 



(5) 



where Mg, Rs are the mass and equatorial radius of 
Saturn, mj,mE the masses of Janus and Epimetheus, 



Table 1: Fundamental frequencies characteristic of the Saturn- 
Janus-Epimethcus system. These frequencies are derived from a 
400 years long numerical integration of the three-body problem in- 
cluding Saturn's oblateness. 





Freq. (rad/day) 


Per. 


n 


9.045924658 


0.69459 (days) 


v 


-0.002147139 


8.01179 (yrs) 


9j 


0.034948139 


0.49223 (yrs) 


9e 


0.034952723 


0.49216 (yrs) 


sj 


-0.034811959 


0.49415 (yrs) 


Se 


-0.034814589 


0.49412 (yrs) 



and a the barycenter of semi-major axis Janus and 
Epimetheus (see Formula (O). The elliptic elements 
(a, e, /, A, TZ7, fij^ are the elements of the unique ellipse 
tangent to the planetocentric velocity at the planetocen- 
tric location of the satellite, assuming that the third Ke- 
pler law is given by equation ([S]). As shown by Figure [TJ 
this definition of elliptic elements removes t he main or- 
bital oscillations from the elliptical elements ( Greenber3 . 
19811 ). In addition, it is well-known that changing the 
value of fi shifts the mean value of the semi-major axis by 
a quantity of the order of J2(i?s/a)^, which leads to the 
translation of about 600 km clearly visible in Fig. m 



1996[) . 3.2 Analytical expression of the elliptic 
elements' variations 



In this section, we detail the quasi-periodic expansion of 
the elliptic elements for both satellites that will be useful 
for the rotation study. 

According to ^ c lassical theor ies (i.e 

Dermott and Murravlll98ll lYoder et al.lll983 or Namouni 
19991 for more recent developments), the variations of 
the mean longitudes and semi-major axes of the co- 
orbital satellites are accurately approximated by the 
expressions: 



Aj ^ \o + nt + CeK mod(27r) 
Xe ~ Xq + ht — C.jXr mod(27r) 



aj ! 
Qe 



; d - Qor 



(6) 



(7) 



where O = rnj/i'mj +™e) and Ce = 1 ~ O- The vari- 
ables ttr and Xr represent the relative semi-major axis and 
mean longitudes of the satellites, that is: = aj ~ oe 



^The index "J" or "E" is added to specify that the elements are 
related to Janus or Epimetheus if necessary. 
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Figure 1: Variation of the semi-major axis of Janus and Epimetheus: 
the upper plot shows the temporal variation of the semi-major 
axes of the two satellites (solid lino for Janus and dotted line for 
Epimetheus), when J2 is included in the third Kepler's law. The 
semi-major axes oscillate around the mean value a = 151436.9 km. 
When J2 is not considered (bottom plot), large short-period oscilla- 
tions are superimposed to the above-mentioned 8-year signal. The 
signals have a period of 277/1/ 8 years. In the latter case, the mean 
value of a is equal to 152024.4 km, which is about 600 km greater 
than when J2 is taken into account. 



and Xr = Xj — Xe- These relations clearly reflect the 
symmetries between the orbits of Janus and Epimetheus. 
Formula ([7]) implies that the barycenter of the semi-major 
axes, namely 0«J + CsdE, is almost constanlH. Numer- 
ical simulation shows that the relative variation of this 
quantity is smaller than 2 x 10""^ if Saturn's oblateness 
is included in the definition of the elliptic elements, while 
the variation is about 200 times more without J2. Simi- 
larly, formula ([6]) implies that the barycenter of the mean 
longitudes Cj^J + CeXe increases almost constantly with 



^Once averaged the Hamiltonian of the three-body problem on 
the mean longitude of the satellites, the quantity Cjy^'V + Ce\/o,e 
becomes an integral of the motion. The relative semi-major axis ar 
remaining always very small with respect to aj and aE, the stated 
property holds. 



time. In other words, {Cj^J + Ce^e) ~ n. It turns 
out that the main variations of a and A are given by the 
relative motion in coordinates (a^., A^). 

Neglecting the terms of powers greater than two in ec- 
centricities (which are very small for these satellites) , the 
relative motion of the satellites satisfies the differential 
systerrd: 



Xr 



2ena (l - (2 - 2 cos Ar)-3/2) sin A^ 

3fi ttr m J + TTlE 

— 7r~j with £ = — — '■ 

2 a AIs + m ; + niE 



(8) 



The solutions of these equations can be expanded in a 
Fourier series as: 



ar(t) — a 



a'j^^ cos i^piyt + 



(9) 



with = 



Xr{t) 



with 



p>i 




;('■) 



sin ( ph't 



4' 



(10) 



where v is the frequency of the relative motion, that is 
the second fundamental frequency defined in Section [Ql 
It corresponds to a period of eight years for Janus and 
Epimetheus. The vanishing of the even coefficients in the 
series ([9]) and (fTO|) is due to the symmetries of the system 
([8]). Indeed, the invariance of the differential system by 
the transformation z 1 — > —z where z = {ar,Xr), leads 
to the relations z{t + T/2) = ~z{t), where T = 271/1/ is 
the period of the solution. Applying this relation to the 



(r) 
2p 



-'2p 



0. 



Fourier expansion of z we obtain a 

While the differential system ([5]) is integrable, it is hard 
to get an analytical expansion of its solutions in the form 

(r) 

of a Fourier series. Consequently, the coefficients aj, , 

(r) (r) 

fip and the phases ipp have been inferred from the so- 
lutions of the numerical integration. A truncated expres- 
sion of this expansion is given in Table [H The compari- 
son between the third and fourth columns, which display 
the coefficients Up and /3p , respectively, emphasizes the 
different decreasing speed of theses sequences. While the 
coefficients of the relative semi-majors axis seem to de- 

(r) 

crease slowly, the sequence (3p converges more rapidly. 
Indeed, Equation (|8]) imposes the coefficients fip to be 
proportional to p~ a), . 

In order to illustrate the convergence of the series ^ 
and (fTU]) towards the solution of the equation of relative 



•^Equivalent formula t ions c an be found in ISalo and Yoded l|l98Sh 
or lRenner and Sicardvl ||2004| ') 
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motion (|S]) , we consider different approximations of tliese 
series for wfiicli the N first terms are summed. To tfie 
relative orbit of tlie two co-orbital satellites in the plan 
{ur, Xr) visible in Fig. [2] (bold curve), is superimposed the 
approximated orbit (dashed curves) obtained by varying 
the integer N. The ellipse obtained for = 1 provides a 
very crude approximation of the relative orbit, while the 
approximation generated at = 30 starts matching ac- 
curately the numerical solution. In addition to the central 
symmetry z i — >■ —z mentioned above. Fig. [2] emphasizes 
a second symmetry with respect to the axis of coordi- 
nates. These symmetries impose relationships between 
the phases ipp, whose description is beyond the scope of 
this paper. 



1 term 

2 term 
5 terms 

10 terms 
20 terms 
30 terms 



-60 -40 -20 20 
aj - ag (km) 



40 



60 



Figure 2: Relative orbit of the satellites: the X-axis represents the 
relative semi-major axis ar = aj — in km, while the Y-axis cor- 
responds to the relative mean longitude Ar = Aj — A^; in radians. 
The solid red curve plots the orbit deduced from the numerical sim- 
ulations, while the dotted curves stands for several approximation 
given by formulas ^ and l)10|l. See text for more details. 

Finally, from the periodic representations © and pUj) 
of the quantity z{t) = (urit), Xr{t)) and according to for- 



mulas (j6]) and the trigonometric approximations of 
the semi-major axis and mean longitude of the satellites 
read: 



\t)=a 



l<p<Af 



(11) 



= a 1 + A 



A(^) (<) = Xa+nt+ J2 ^p""^ sin (piyt 



Xo+nt + B 



1<P<N 



N 



Pl^^ = {I - cM'^K 



(12) 



(13) 



where the index x replaces J and E, whether wc consider 
Janus or Epimetheus. The numerical values of the coef- 



ficients ttp^"* and /3p'^' are reported in the fifth and sixth 
columns of Table[2]for Janus and in the seventh and eighth 
columns for Epimetheus. In addition, the ratios of each 
ap^''/ap'^'' and /Jp^'/Z^p'^'' listed in Tableware close to 3.6, 
in agreement with the formulae (|13p. The short-period 
oscillations do not appear in that table because they are 
negligible in comparison to the other parameters. 

In the following section, we use this representation of 
the elliptical elements of the satellites, especially for the 
mean longitudes, to develop an elementary perturbation 
theory describing the rotation of Janus and Epimetheus. 
We consider in the next section that the eccentricities and 
inclinations are constant, as underlined in section [XT] and 
discussed in Section 14.21 In the same way, the preces- 
sion of the periccntcrs and nodes will be approximated by 
assuming uniform motion defined as: 



e(^) ~ e, zu{t) = gt + zuq 

i{t) = /, n{t) = st + no 



(14) 
(15) 



4 Physical librations for Janus 
and Epimetheus 

4.1 Perturbative analysis 
4.1.1 Dynamical equations 

The equation governing the physical libration is inferred 
from the angular momentum balance equation projected 
onto the equatorial plane of the body 

e + l^^sin2i0^.)^, (16) 
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Table 2: Periodic approximation of the semi-major axes of Janus and Epimethcus. The relative quantities a,- and \r according to Equations 
l(9t and 1)1011 arc represented in the third and fourth columns. The fifth and sixth columns contain the coefficients of a^''^ and A^"') appearing 
in Equations l lllll and II12II . while the seventh and eighth columns correspond to the same quantities for Epimetheus. The frequencies listed 
in the second column are equal to pv where p is in the first column and v is given in Table [T] The phases in the ninth column apply directly 
to Equation l(9]l and l llOII , and can be adapted to Equation l lllll and I I12II by means of Equation I I13II . The origin of time in formulas is 
equal to 1949-Dec-28 00:00:00.0000 (JD 2433278.5) from the ephemeris Horizons. 
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where all variables have been defined in Section 2. Since 
the angle 9 — v remains always small, the linearization of 
the equation (|16p is a valid approximation to the rotation. 
Therefore, in this section, we will consider the linear time- 
dependent equation: 

a\^/. X ,9 GMs(B-A) , , 

- ef-w =0, with = 3 — ^ ' 17 
r J (J 

This equation is not integrable because the quantities ct, 
ajr and v are implicit functions of time. However, by 
using the expansions in quasi-periodic functions of time 
presented in section |3] for the orbital elements, and an el- 
ementary perturbation theory, we can obtain the approx- 
imated solution of this equation. Let us first introduce 
the physical libration 7. This angle is defined as the os- 
cillation of B around the uniform motion (n — s)t, where 
n — s is the main frequency of the draconic mean longi- 
tude A — ri = £+uj. Therefore the physical libration reads: 
7 = — (n — s)t — £q — loq, the angles £q and loq being the 
initial values of the mean anomaly and of the argument of 
the pericenter of the considered satellitcQ. Consequently, 
using the relations p2|) . p4|) . (|T5l) . the definition of 7 and 
V = f + io for the true anomaly /, the angle 6 — v also 
reads: 

6 — v = j — v + {h~ s)t + ^0 + ^0 

= -/-[£ + uj-in- s)t -£o- ujo] -[f-e] (18) 
= 7 — Bn — 2esin^ 

where terms of order 2 and greater in eccentricity have 
been neglected. At this point, it is convenient to use the 
function y ~ ■y ~ Bn. Indeed, although the amplitudes of 
the terms contained in Bn are very large, their frequen- 
cies are small and therefore the acceleration generated by 
Bn, which is of order j/^, is negligible with respect to a^. 
Under these approximations. Equation (|17p becomes: 



y 



(19) 



This equation has the same form as in the Keplerian case 
(see Equation (jSj), but it is not periodic anymore but 
quasiperiodic on the time because, according to formu- 
las ((HI) and ((HI), £ reads: 



£^X-w ^£o + ht + B 



with 



so 



i£ iin int 

e = e 'e 



n e 

l<q<N 



l<<j<7V kez 



(20) 



(21) 



^The definition of 7 is different from the one used in the Keplerian 
case for in the present situation the angle £ + ui is not proportional 
to the time. 



where the Jk{x) are the Bessel functions (see Ap- 
pendix]^. Let us recall that in the previous expression, 
the coefficients /3g vanish when q is even. Applying usual 
properties of the Bessel functions that are recalled in Ap- 
pendix [XI wc deduce that, for q odd, the ratio between 
the coefficients Jk{Pq) computed for Epimetheus and for 
Janus is well approximated by: 



Jfe(/3f^) 



(22) 



For this reason, the coefficients of the expansion (PT|) de- 
crease with increasing k much more rapidly for Janus than 
for Epimetheus. Consequently, the number of terms nec- 
essary to approximate e*^ to a given accuracy using a trun- 
cated expression of ([21]) is different for the two moons (see 
Section l4.2.2|) . In order to simplify the following develop- 
ment, we present the series containing only the first term 
A'^ = 1 , which corresponds to approximate the relative or- 
bit of the two bodies by the green ellipse in figure [2] This 
makes a crude simplification regarding the orbital motion 
of the moons but the generalization to N higher than 1 is 
then straightforward. Then we have: 



(23) 



fcGZ 



It turns out that, under these approximations and assum- 
ing as in ([T4| that the eccentricity is constant and denoted 
e: 

e sin £ = eS with 

S= Jo(/3i)sin(nt + 4) + 

^JpiPijli^i-i^iin+pv^t+pipi +£o)+ (24) 
p>i 

(-l)P sin((n - pv)t ~ pipi + £„)] 

Let us mention that in the expression (j23p the index of 
the summation k belongs to Z, while in ([24| the summa- 
tion is restricted to positive integers, p. Using ([231) and 
expanding (a/r)^ at first order in eccentricity, we also get 
the expression of: 



- ) = 1 + 3eC 



(25) 



The expression of C is the same as iS, where sine functions 
are replaced by cosine functions. The last term that we 
have to expand is = <T^{a/a)^. By ([TT|) with A^ = 1, 
wc have a'^ w a'^ (1 -f 3^i) and consequently: 



(1 - 3A) with 0-2=3 



GmB-A 



(26) 
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By substitution of the relations and in by the expression: 

the expression (jl9p this equation becomes: 



y+a^{l-3Ai){l+3eC)y = 2ea^{l-3Ai){l+3eC)S (27) 

If we spht y into a sum of terms of decreasing magnitude 
as y = yo + yi + ■ ■ ■ , we obtain the fohowing system of 
equations: 



(28) 



yi + CT^yi = 3ct^ {Ai - eC) 2/0 + 6eCT^ (eC - yti) S (29) 



The solutions of these equations arc trigonometric sc- 
ries, whose frequencies are linear combinations with in- 
teger coefficients of fundamental frequencies of the satel- 
lites (n, v, g, s) and of the frequency of the free libration 
a. The solution independent of the free frequency a is 
usually called "forced solution" and it is a quasi-periodic 
series of frequencies (n, v,g, s). 

When weak dissipation is introduced in the system, al- 
most tra jectories converge towa r ds a q uasiperiodic attrac- 
tor (see ICelletti and Chierchial (|2008l )). This quasiperi- 
odic trajectory is very close to the forced solution de- 
scribed in the conservative system. Indeed, the funda- 
mental frequencies of these two solutions are the same, 
only the amplitudes and phases are slightly different (see 
section [S^ . Consequently, it is relevant to focus on the 
forced solution in the conservative problem. 



4.1.2 Forced librations 

To begin with, let us associate to a quasi-periodic function 
/ the function / such that: 

p 

then f{t) = ^ -.j _ ^ fp sin{vpt + (j>p) 



With these notations, the general solution of ([28]) reads 



yo{t) = hsu-i{at + V) + 



(30) 



where h and ^ are arbitrary constants. In this section, 
we focus on the forced solution, so we put h = 0. As a 
consequence, It is easy to verify in ([28]) that the contri- 
butions of eyoC and e^CS, denoted respectively ey^C and 
e^CS^, is a second-order in eccentricity, and that Aiyo and 
eA\S are of order eai. As, for Janus and Epimetheus, the 
coefficient ai is lower than e, then the term yi can be ne- 
glected. Finally, the forced libration can be approximated 



l<q<N 



+ 



2eaVo(/3i) . 



sin {fit + io)+ 



2ea2^Jp(/3i) 



p>i 



sin((n + pi^)t + pfi + £q) 
— {n + pv)"^ 



(31) 



^ -^^p sin((n - piy)t - pipi + io) 
— [h — pv)^ 



Let us recall that, in order to give a simple expression of 
the short-period libration, we put iV = 1 in the expression 
((2T|) . For this reason, only /3i appears in the short-period 
part of pip . In contrast, in the long-period component 
of 7 (first term of Formula (|3T|) ). N is arbitrary. The 
librational angle 7 is split in two types of terms exhibit- 
ing different behaviors. The first type corresponds to the 
27r/i^-periodic terms that depend only on the coefficients 
/3p, i.e., on the mean longitudes of the satellites. For these 
long-period terms, the dynamical figure has no influence. 
The second type includes terms that vary rapidly (quasi- 
periodic with short frequencies h±pv) and depend on the 
triaxiality of the body {B — A)/C through the libration 
proper frequency a. 

The amplitudes of the rapidly oscillating terms depend 
on the magnitude of the forcing, 2eJp{pi) and on the prox- 
imity of the forcing frequency h zL pv with the libration 
proper frequency a. In the case of Janus, the proper fre- 
quency is 4.96 rad/day, which is far from the resonance, 
whereas for Epimetheus the proper frequency is equal to 
8.52 r ad /da ys, and its influence on the amplitude is sub- 
stantial (see iTiscareno et aL ( 20091 )). 

By contrast to the Keplerian case, the swap results in 
the amplitude of the term associated to the frequency 
n = n — 5 to be proportional to Jo(/3i). This term is of 
the order of 1 for Janus, but it is significant in the case of 
Epimetheus as close to 0.22. Therefore, for both satellites, 
the rotation significantly departs from the Keplerian case. 



4.1.3 Proper libration 

In this section we investigate the proper libration of the 
moons (also called free libration) and we especially focus 
on the influence of a small divisor on the solution. By dis- 
sipative effect the proper libration is expected to be small 
and the damping time scale is short. However excitation 
mechanism might exist like recent impact for example. 
To focus on proper libration, let us remove the external 
forcing by imposing e = in (|28l) and (|29l) . Then, the 
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solution (j30p reads yo{t) — hsm{at + ip), where the am- 
plitude h is small but different from zero. Consequently, 
by substitution of yo in the equation (j^ . we get: 



3a^aih f sm{{u + a)t + ipi + ip) 
sin((i^ — a)t + (^1 — -0) 



3 aa\ 
2~' 



(32) 



Contrarily to the case of the forced libration, the ampli- 
tude of the term Ai, which is equal to ai (see Tabled, 
plays a major role here. Indeed, in the present situation, 
that term is multiplied by the factor ct/i^, which is about 
2200 for Janus and 3700 for Epimetheus. Thus it gener- 
ates a second-order solution yi , whose size is comparable 
to the solution of order one y^. Then, using the values of a 
given in Table[3]for Janus and Epimetheus, the librational 
responses for the two satellites are: 

yj = h (sin((Ti + V') + 0.3sin(i^f + ipi)cos{dt + i/")) (33) 
yE ~ h (sin(tTi + ip) + 1.87sin(i^i + ipi)cos{dt + ijj)) 

(34) 

The proper librations are combinations of a sine term with 
a constant amplitude and a cosine term with an amplitude 
varying at the swap frequency. For Janus, the main term 
is the sine component, whereas for Epimetheus, it is the 
cosine component. 

4.2 Numerical study and quasiperiodic 
representation of the libration 

4.2.1 Quasiperiodic decomposition 

We numerically integrate the orbital and rotational dy- 
namical equations (|16|) with the triaxialities listed in Ta- 
bic [3] in order to determine the accuracy of analytical so- 
lution. To clearly separate in the frequency analysis the 
frequencies (n ±pj/), which are quite close, we have to 
integrate the trajectories over a long time-span of about 
400 years in the future and 400 years in the past. We 
also integrate during 8000 years to study the very slow 
variations of the eccentricities of the two satellites. This 
point will be discussed below. 

Wc focus on the forced libration but the initial condi- 
tions of such trajectory are not known. Hence, we use an 
iterative method based on the frequency analysis to con- 
verge towards this trajectory by re r novin g the free libra- 
tion amplitude (see Couetdic et al. ( 2009( ) Section 4.5). 

As predicted by the theoretical approach stated in sec- 
tion I4.1.2[ the libration angle 7 = 9 — fit — £0 — u>o 




500 1000 1500 2000 2500 3000 3500 4000 

t (years) 

Figure 3: Secular variations of the eccentricities of the two sateUites 
(sohd curve for Janus, dashed curve for Epimetheus) due to their 
mutual gravitational interactions. The grey rectangle represents the 
time interval used for the numerical simulation of the rotation. 



can be naturally written as the sum of two components: 
7 = 7; -f 7s , where 7; contains the long-period harmonics 
and 7s the short-period ones. The frequencies of these 
harmonics take the form: h ± piy. The long-period oscil- 
lations whose amplitude is about 36° for Janus and 130° 
for Epimetheus, are very well described by the first sum- 
mation in the expression pip, that is: 



7i(*) = X! /3gSin(gi/t 

g>l, q odd 



where the numerical values of the coefficients /3q appear 
in the sixth column (for Janus) and in the eighth one (for 
Epimetheus) of Table [2] In other words, the Fourier ex- 
pansion of 7/ is similar to the one governing the long-time 
oscillations of the mean longitude of the corresponding 
satellite. 

In order to estimate the accuracy of the coefficients of 
the short-period terms obtained analytically in Section 
14.1.21 we will focus on the Fourier decomposition of func- 
tion 7s deduced from the numerical simulation. 

According to the equation ([T9|) . the amplitude of the 
short-period component of the rotational libration 7s is 
proportional to the eccentricity of the satellite. As this 
quantity suffers from large very long-period variation (pe- 
riod of about 2760 years), the amplitude of 7s is slowly 
time dependent. The evolution of the eccentricities are 
represented in Fig. [3l During the 400-year numerical 
integration, the variations of the eccentricities are sig- 
nificant (according to Fig. [3J e,/ increases by ~ 20%, 
while ce decreases by more than 40%). Consequently, 
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Table 3: Physical properties of Janus and Epimetheus. (1) Tiscareno et al. (2009). 



Janus Epimetheus 

Mean Radius (i)(km) 89.5 ± 1.4 58.1 ± 1.8 

Long axis(i) (km) 101.5 ± 1.9 64.9 ± 2.0 

Intermediate axis^^) (km) 92.5 ± 1.2 57.0 ± 3.7 

Small axis(i) (km) 76.3 ± 1.2 53.1 ± 0.7 

Triaxiality(i) (B-y4)/C 0.100 ± 0.012 0.296t5!-^27 

a (rad/day) 4.95508166 8.52504483 
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Figure 4: Relative error €p of the Fourier approximations of 7 
against the number of terms included in the decomposition. The 
" + " correspond to Janus and the " X " to Epimetheus. 



Figure 5: Relative error of the Fourier approximation of 
7 in the interval [0 : 200] years = [JD 2433278.5 : 
JD 2433278.5] for Epimetheus. The instantaneous relative error 
(7s (i) — 7s (t)) /Max|7s(t)| is plotted in red for p = 13 and in green 
for p = 57. 



the amplitudes of the short-period librations, which are 
proportional to the eccentricities, are modulated with 
the same ratio. As a consequence, we study the quan- 
tity 7s(t) = 7i5(i)e(0)/e(i). The multiplicative factor 
e(0)/e(t), ratio of the initial eccentricity by the eccen- 
tricity at the time t, imposes to the amplitude of 7s(i) to 
be almost constant during the integration time. During 



the 200 first years of the simulation, the ratio e{t)/e{0) 
can be accurately fitted by the quadratic polynomial P, 
equal to P jit) = 1 + 1.9839 x lO^^i - 4.1467 x 10" ^^^^ foj. 
Janus, and PE(t) 1 - 2.3486 x lO^^t - 1.7633 x lO^^^^^ 
for Epimetheus. The time is counted in Julian days from 
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1949-DCC-28 00:00:00.0000 (JD 2433278.5). These poly- 
nomials approximate the considered ratio with a relative 
accuracy of 0.32% for Janus and 0.45% for Epimetheus. 
Consequently, the value 7^ {t) of the short-period compo- 
nent libration at an arbitrary time t inside the 200-years 
considered interval can easily be deduced from 7s(i) using 
the relation js{t) = P{t)js{t). 

The main terms of the Fourier decompositions of 7^ are 
displayed in Table|3]for Janus and Table[5]for Epimetheus. 
The solutions are given in the form: 



7(i) = sin{fpt + ipp) 

p 

with fp — jph + kpV 



(35) 



The first and second columns contain the amplitudes of 7^ 
from the numerical simulation and the analytical expres- 
sion ()3ip . The comparison between these two columns 
will be discussed in the next section. The third column 
contains the frequencies /p, while the integers jp and kp 
are displayed in the fourth and fifth columns. The last 
column presents the phases il^p deduced from the numer- 
ical simulation. 

In the case of Janus (Table |4]), the term at the orbital 
period is dominant, with an amplitude of 0.0052 radians, 
that is 0.3°. By adding the following terms the amplitude 
is slightly modified. We deduce that, from the polyno- 
mial interpolation of ej given above, the amplitude of 7s 
increases for 0.34° in 1950 to 0.36 ° in 2010, which is ver y 
close to the value of 0.33° given in lTiscareno et aL ( 20091 ). 

For Epimetheus (Table [5]), no dominant term appears 
clearly and the decrease of the coefficients in slower than 
for Janus. The amplitude of the short-period oscillations 
given by the sum of this 13 terms deceases from 8.77° in 
1950 down to 8.35° in 2010, which is here ag ain, compa- 
rable the 8.9° given in lTiscareno et all (|2009f ) (Table 7). 



Let us now briefly discuss the accuracy of the quasiperi- 
odic approximations of the numerical solutions. Due to 
windo wing used in the frequency analysis (see iLaskar 
([20051)), the method docs not provide a uniform approxi- 
mation of studied signal. The accuracy is generally lower 
on the margin of the considered interval of time. Con- 
sequently, the time span of the analysis has been chosen 
such that the best accuracy is obtained in the interval 
/ = [0 : 200] years. On the interval /, the accuracy is 
measured with the help of the relative error €p defined as: 



Max|7,(t)-7f(t)| 
Max|7,(t)| 



(36) 



In his expression, 7I' is the quasiperiodic approximation 
of 7<, containing the p dominant terms of the decompo- 



sition. As mentioned above, the value of Max|7s(t)| is 
about 0.34° for Janus and 8.77° for Epimetheus. 

The number of terms given in Tables 2] and [5] is such 
that the relative error Cp is better than 10% for both satel- 
lites (Janus 8.7%, Epimetheus 9.3%), which gives an ab- 
solute error of about 0.8° for Epimetheus and less that 
0.03° for Janus (The number of selected terms is 13 for 
Epimetheus and 7 for Janus). The increase in the ac- 
curacy with respect to the number of terms included in 
the quasiperiodic approximation is presented in Fig. m 
where logiQ^Cp) is plotted against p. As it was already 
mentioned, the convergence of the Fourier approximation 
is much rapid for Janus than for Epimetheus. This figure 
also shows that 45 terms for Janus and 57 for Epimetheus 
are necessary to reduce the relative error down to 1.2%. 

Finally, Fig. [5] shows the instantaneous relative error 
{ls{t) - T-f (i)) /Max|7s(t)| for p = Yd (red dots) and p = 
57 (green dots) in the case of Epimetheus. The results are 
similar for Janus. For p = 13, the discrepancy is mainly 
due to an over estimate of the influence of the orbital swap 
on the rotation (one peak every four years). The addition 
of the harmonics of frequency v ± pfi for high values of p 
tends to erase this effect. 



4.2.2 Comparison of the analytical and numeri- 
cal QP representations 

The analytical and numerical short-period component of 
the forced solution are listed in Table 2] for Janus. The 
accuracy of the short-period terms is a function of the 
considered harmonic, whose frequency is n ± pu, where p 
is an arbitrary positive integer. For p — the discrep- 
ancy is about 3% and increases for increasing p, e.g. 14% 
for p = 2 and to 85% for p = 2>. This lack of accuracy 
can be ascribed to neglecting the terms Pq for q strictly 
greater than one in the expansion (|2ip . and consequently 
in the analytical solutions of the rotation pip . Indeed, as 
it can be shown by a straightforward calculation, the am- 
plitudes of the terms that have been neglected are given 
by 2J3(/3i)Ji(/33) ~ 2 X lO'^ for p = 0, J2{(ii)Ji{(is) « 
4 X 10-3 forp= 1 and Jj (/?i) Jj (/?3) w 2 x 10"! forp = 2. 
These numerical values, which are deduced from Table [6l 
are in good agreement with the level of accuracy men- 
tioned above. 

The librational behavior of Epimetheus is reported in 
Table [5) Contrarily to Janus, the amplitudes of the 
short-period terms at n ± v and n ± 2iy are greater 
than the terms at n. This is due to the fact that 
Jiif^i) > J2{Pi) > Jo{(3i) for Epimetheus, while the rela- 
tion Jo(/3i) > Ji(/3i) > J2{Pi) holds for Janus (see Tables 
[S] and [T]) . From comparing the first two columns of Table 
[5] we find that the accuracy on the short-period terms ob- 
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Table 4: Frequency analysis in the short-period librational motion 
of Janus. This table displays the 7 first terms of the quasipcri- 
odic decomposition of 7(4) = e(0)/e(t)7s (t). See the text for more 
details. 



Amp num 


Amp ana 


Freq 


nj 


u 


phase num 


(rad) 


(rad) 


(rad/days) 






(rad) 


0.0052021 


0.0053573 


9.0109765 


1 





4.705060 


0.0014761 


0.0015412 


9.0131237 


1 


-1 


3.151328 


0.0014756 


0.0015433 


9.0088293 


1 


1 


3.117592 


0.0002567 


0.0002159 


9.0152709 


1 


-2 


1.596407 


0.0002521 


0.0002165 


9.0066821 


1 


2 


1.532845 


0.0001427 


0.0000201 


9.0045349 


1 


3 


3.087729 


0.0001412 


0.0000200 


9.0174181 


1 


-3 


3.183252 



when N, in the series S, tends towards the infinity. To 
solve (28) we have expanded y in Fourier series, here we 
use the fact that two time scales coexist in the problem in 
order to solve the equation (37). The frequency 1/ being 
very small with respect to n, it is possible to expend the 
solution in Taylor series of the small parameter e = v/n. 
First, let us remember that, according to (PH]). the mean 
anomaly £ can be written: £ = nt + B{vt) were B contains 
the long-time variations associated to the orbital swapH. 

If we now introduce the three angles {0i,02, ^3), defined 
by Oi = nt, 02 = vt and 0^ = d't, we can see any solutions 
of ([57)) . as a quasiperiodic function of the three previously 
defined angle. More precisely, this equation being linear, 
its solution reads: 



tained for Epimetheus is worse than in the case of Janus. 
It reaches about 9% for p = 1 and increases to 28% for 
p = 4 and even more than for 80% for p > 5. As for 
Janus, the accuracy of the analytical solution would be 
strongly increased if the terms related to /^a and possibly 
to 13^ were taken into account. 

Table 5: Frequency analysis in the short-period librational motion of 
Epimetheus. This table displays the 13 first terms of the quasiperi- 
odic decomposition of "/{t} = e{0) / e{t)^s{t) ■ See the text for more 
details. 



4.3 



Adiabatic-like solution of the equa- 
tion of the hbration in longitude 



We have seen in section 14.1. 1[ neglecting the quadratic 
terms in eccentricity, that the short-period component of 
7 can be approximated by Eq. ([28]) , which is the equation 



2ea sinl 



(37) 



y( 



h) = y{0iM + v{Oi) 



(38) 



y being periodic and y quasiperiodic. The forced solu- 
tion is consequently y. Therefore, in order to obtain this 
particular solution, it is enough to restrict our analysis to 
the set of the quasiperiodic functions of the two variables 
(01,02)- Consequently, if y is an element of this set, its 
time derivative reads: 

dy _ dy dy 



dt "a^i '^'''802 



(39) 



It turns out that the forced solution of the equation (|57)> 
is also solution of the linear partial differential equation: 



Amp num 


Amp ana 


Freq 


flE 


V 


phase num 


(rad) 


(rad) 


(rad/days) 






(rad) 


0.0816497 


0.0892555 


9.0088248 


1 


1 


1.356172 


0.0813811 


0.0884486 


9.0131192 


1 


-1 


1.389845 


0.0618972 


0.0536462 


9.0152664 


1 


-2 


2.977836 


0.0617296 


0.0546295 


9.0066776 


1 


2 


2.908991 


0.0405482 


0.0349633 


9.0109720 


1 





6.085400 


0.0172653 


0.0194522 


9.0174136 


1 


-3 


4.566668 


0.0168758 


0.0199895 


9.0045305 


1 


3 


4.462410 


0.0073675 


0.0052782 


9.0023833 


1 


4 


2.877236 


0.0072106 


0.0050899 


9.0195608 


1 


-4 


3.006457 


0.0052737 


0.0010462 


9.0217080 


1 


-5 


4.596183 


0.0051221 


0.0010948 


9.0002361 


1 


5 


4.429255 


0.0024037 


0.0001873 


8.9980889 


1 


6 


2.845404 


0.0023006 


0.0001774 


9.0238552 


1 


-6 


3.038194 



-2^ 

" 801 



2hv 



d'^y 
d0id02 



2^ 

' 901 



a^y 



(40) 



2a^esm{0i+B{02)) 



If we now expand the solution y in Taylor series with 
respect to e, that is: y = yo + eyi + • • • + e^yp + • ■ • , the 
functions yp satisfy the sequence of differential equations: 



2/0 
2/1 
Vp 



.2^^2/0 



" 801 

-2^^2/1 



0-^2/0 



801 
801 



2" + 2/1 



25-2esin(6'i +6(6*2)) (41) 
9^2/0 



c^^2/p 



80i802 

80i802 
2" 



(42) 
(43) 



n 



801 



Vp > 2 



Solving iteratively these equations, it is easy to show 
that the forced solution of the libation equation reads: 

2/(0 = ^ [(1 + £U{vt, e)) sin{nt + B{vt)) 



(44) 



e^Vivt^e) cos {nt + B{vt))' 



''The subscript N has been removed from B to indicate that the 
summation can be infinite 
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where U and V are periodic functions in that can be 
expanded in power series of e as: 



U 



p>Q 



and V 



p>0 



(45) 



The coefficients Up and Vp, which depend on B and on 
it first p derivatives, are deduced form (|4T|) to (|43)) by 
induction. The approximation of zero order (when v is 
neglected with respect to h): 



yo{t) = sin {ht + B{vt)) _ 



■ sm< 



(46) 



is the solution that we intuitively get freezing the long- 
term temp oral variations. This s olution is close to the 
one used in lTiscareno et ahl ( 2009f ) for the libration, that 
is: Asnv(n{t) t + (p), where both amplitude A and phase 
(j) have been fitted to the observations. In this expression, 
the instantaneous mean motion n{t) reflects, like B{vt) in 
formula pS)) . the long-term variations of the "instanta- 
neous" mean motion due to the orbital swap. 

Even if the solution yo given by (|46|) resembles the li- 
bration in the Keplerian case (section [5]) , their behaviors 
are very different. Indeed, in the Keplerian case the mean 
anomaly l increases linearly with the time, while in the 
more realistic case that we are considering, a large 8-years 
periodic motion in superimposes to this linear evolution 
(about 36° for Janus and more than 130° for Epimetheus). 
Forget this long-term variation could lead to an error on 
the amplitude of the libration reaching at worst 50% for 
Janus and 100% for Epimetheus. 

In order to evaluate the accuracy of the "frozen" solu- 
tion j/o given by (j46|) . we have to estimate, at least, the 
size of the term Ui which appears in the expression of yi. 
Substituting yp and solving the differential equation (|42|) , 
we have: 



v{t) 



Tl \ 

1 + 2e—^ T^B'{vt) X 

0-2-^2 ^ 'J (47) 



sin {fit + B{vt)) + 0{e^) 



where B' is the first derivative of B. 

Deducing from numerical simulations that the upper 
bounds of \B'\ are respectively 1.5 for Epimetheus and 
0.45 for Janus, it turns out that the addition of yi in the 
forced solution, modified the amplitudes of the "frozen" 
solutions of about 0.64% for Epimetheus and 0.03% for 
Janus, which is far lower than the accuracy of the obser- 
vations. 



5 Discussion 

5.1 Higher harmonics 

Tiscareno et"!!! (|2009f ) found in the shape fitting residu- 



als an unexplained offset in the direction of the longest 
figure axis of the moons, 5.2° in the case of Janus. For 
Epimetheus, the detected offset is within the error bars. 
These authors suggested that the departure observed for 
Janus is due to large lateral density anomalies. We inves- 
tigate such a hypothesis by assuming that the satellites 
shapes depart from triaxial, hydrostatic shapes, due to 
mass anomalies expressed at the thi rd degree of sp heri- 
cal harmonics (as for the Moon e.g. lEckhardt 1981 ). In 
this case, in addition to periodic terms and small shifts 
in the proper frequency, a constant term appears in the 
dynamical Equation leading to a constant offset of the 
form 

7hs3 = 4(-15^33 + 0.5531) (48) 



where C = (^j (f ) = 0.308 rad/days^ for Janus and 
0.198 rad/days^ for Epimetheus, where 533 and 531 are 
spherical harmonics of order 3 and degree 3 and 1, re- 
spectively. In the expression p8)) . we have assumed that 
the spherical harmonics C13 and C33 are negligible with 
respect to €22- Thus, the offset can be expressed as 



7hs3 



for Janus and 



7lfs3 



0.72°(-15533 + 0.5531) 



0.16°(-15533 + 0.5531) 



(49) 



(50) 



for Epimetheus. 

As a consequence, an offset of 5.2° requires that the 
combination (— I5533-I-O.5531) is about 7 for Janus, which 
seems very large even for rubble-piles. In this case, the 
shape would be far from an ellipsoid. Nevertheless, the 
contribution of density anomalies cannot be completely 
ruled out, and more complex shape models remain to be 
developed in order to better assess the infiuence of non- 
hydrostatic anomalies on rotation. 

5.2 Tidal dissipation 

A second pos si ble o rigin of offset determined by 



Tiscareno et al 



(|2009t ) might be related to the tidal 
torque. Indeed, Saturn raises a tidal bulge on each moon 
that is shifted from the planet-satellite direction due to 
the inelastic response of the moon material. Saturn ex- 
erts a gravitational torque on this tidal bulge and the 
body responds by displacing its permanent bulge so that 
it cancels the average saturnian torque acting on the tidal 
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bulge. Such a displacement has been measured for the 
Moon (Williams et al. 2001) and estimated for Enceladus 
(Rambaux et al. 2010). First, we analytically introduced 
the main tidal deformation as a supplementary term in 
equation (|37|) . describing the adiabatic behavior of the 
moons, of the form 



(51) 



where R is the satellite's radius, fc2 the tidal Love num- 



ber, and = {f) 



UiUj. The Ui are the direction cosines 



between Saturn and the satellite. The symbol star means 
the position of Saturn at the constant time delay St re- 
sulting from the dissipation. By taking advantage of the 
spin-orbit synchronous resonance the quantities ui and 
ul are of the order of the unity and U2 and ^2 are small. 
Then, we develop U2 in Taylor series with respect to St. 
Finally, the deformation included in the adiabatic equa- 
tion Eq. (|4T|) reads 



(52) 



ij + 2Xy + a^y = 2ea^ sm£ + ienXcosi 



with the dissipative rate 



2A 



C Gm 



-St 



(53) 



where m is the mass of the satellite. We assume that the 
delay is constant and equal to: St = n~^Q~^. Following 
the perturbative analysis performed in section I4.1.1[ we 
solve this equation by written the particular solution in 
the form 

y ^ ys sin £ + yc cos £ 

where ys is the amplitude of the in-phase term and yc is 
the amplitude of the out-of-phase term raises by dissipa- 
tion. The last term induces a small displacement each 
time that £ = mod(27r), of 



^1^ 



Vc = 



-S2"" 



with 



and 



6e fh^R^ 



C V Gm 



9 /n'^R^ 



C2 V Gm 



(54) 



(55) 



(56) 



The coefficients Si and S2 are equal to —0.011 and 0.3205 
radians for Janus and —0.603 and 11.863 radians for 
Epimetheus. j/c depends on the rheology of the satellite 



through k2/Q. This ratio depends on the internal struc- 
ture and therefore on the origin of the bodies. Charnoz 
et al. (2010) have suggested that many small satellites of 
Saturn, and especially Janus and Epimetheus, come from 
the accretion of ring material in the form of lumps that 
separate from the rings. We expect this accretion scenario 
to yield homogeneous bodies whose composition is mostly 
water ice. For the two small moons the Love number ^2 
may be computed from the following relationship (e.g., 
McDonald 1964): 



1 



21 



2pgR 



(57) 



where p is the mean density for the satellite, R its mean 
radius, g the average surface gravity, and Hr the effec- 
tive shear modulus characteristic of the rubble material. 
We infer from Goldreich and Sari (2009) that the effec- 
tive shear modulus of Janus is about 0.3 GPa and that 
of Epimetheus is 0.2 GPa. Such a low effective modu- 
lus implies an increased value of k2 with respect to that 
expected for a monolith, of the order of 4.8 x lO^"' and 
3.3 X 10"'' for Janus and Epimetheus, respectively. 

With regard to the dissipation factor, its value for ag- 
gregates is poorly constrained and depends on the dis- 
sipative mechanism acting inside the satellites. Then a 
dissipation factor of the order of 10 to 10^ is a possible 
range of values for Janus and Epimetheus. The smallest 
dissipation factor [Q = 10) leading to the largest displace- 
ments of the axis of figure of both satellites and implies 
an yc about 3.02 x 10~^ degrees for Janus and 1.14 x 10~^ 
degrees for Epime theus, which is too sma ll to explain the 
offset observed by iTiscareno et al. I (I2OO9I) . 



5.3 



Influence of the triaxiality on the h- 
brational amphtudes 



The short-period librations are of geophysical interest be- 
cause their amplitudes depend on triaxiality, as shown in 
Eq. (PT|) . We range the triaxialit y inside the error bars 
provided bv lTiscareno et al.l (|2009l) for Janus 0.100±0.012 
and for Epimetheus [0.269 : 0.315]. For Janus, the libra- 
tion amplitude depends linearly on triaxiality, while for 
Epimetheus this dependence is hyperbolic because the 
proper period and the orbital period are close to each 
other (see Fig. [S]). 

The model suggested by [Ch arno i (I2009I) leads to ho- 
mogeneous bodies. Alternatively, Porco et al. (2007) 
suggested that both satellites could have accreted rings 
particles around a core of satellite material with a lower 
porosity. This idea is supported by the fact that accretion 
models can account for the very oblate shapes of the satel- 
lites. That model implies a contrast in density between 
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0.45 




I 0.25 - 




(B-A)/C 



Figure 6: Influence of the triaxiality on the amphtude of the mode 
n — (red curve) where k is equal to J for Janus (a) and E for 
Epimetheus (b), n — gj^ + u (green curve), and n — gi^ — v (blue) 
from Eq. (31). The black bold curves correspond to the amphtude 
deduced from the adiabatic invariant model (Eq. I46p . 



6 Conclusion 

In this paper wc have investigated the hbrational motion 
of the co-orbital satcUites Janus and Epimetheus by us- 
ing three methods: (1) a perturbative technique based 
on quasi-periodic expansions, (2) an adiabatic invariant 
approach by expanding in power series of the small pa- 
rameter V /n, and (3) a numerical integration. With the 
perturbative technique, we have detailed the hbrational 
behavior. For both satellites the solutions are composed 
of long-period librations linked to the orbital swap and 
short-period librations related to the orbital period. We 
found that the amplitudes of the short-period librations 
depend on the magnitude of the forcing and the proxim- 
ity to the resonance, as for the librations analyzed in a 
Keplerian framework, but also on Bessel functions of the 
amplitude of orbital libration of the moons mean longi- 
tudes along their horseshoe orbit. These amplitudes bear 
the signature of the mass distributions in the satellites 
and are crucial to investigate the internal structure signa- 
ture of these objects. On the other hand, the amplitudes 
related to the long-period librations do not contain any 
information on the distribution of mass. The numerical 
integration allows us to assess the accuracy of the pertur- 
bative development. The accuracy of the analytical solu- 
tion is good for Janus but poor for Epimetheus because its 
mass is smaller than Janus and therefore its dynamics is 
more perturbed. In addition, in order to obtain a compact 
analytical solution easy to manipulate, we have developed 
an adiabatic approach yielding directly the amplitude of 
the short periods. 

The analytical approaches have been developed in the 
most general formalism and may be applied for co-orbitals 
like Telesto, Calypso, Helene and Polydeuces. The adia- 
batic approach seems a convenient approach to fit the 
short-period librations to the observations. 
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A Expansion in Bessel functions 

The Bessel functions can be defined as Fourier's coeffi- 
cients of the 27r-periodic function u i — > ^ixsmu -^^j^ere x 
is a real parameter, that is: 

k—~oo 

1 f^^ 

with Jk{x) = — / expi{x sin u — ku)du 
27r Jo 

1 f^^ 

= — / cos(a:sinM — fcu)(iM 
27r Jo 
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These functions satisfy the two following relations that we 
use in Section 14.1.11 : 



Table 7: Numerical values of the coefScients Jq(l3p) in the case of 
Epimetheus 



for allp G 



for allp e N, 



J-p{x) 
Jp{x) = 



{-If 



2Pp. 



,(1 



Jp{x) 
0{x^)) 



In addition, we show in Tables [6] and [7] the values of 
the main coefficients Jq{fip) that are greater than 10^^. 
These tables are useful to evaluate the accuracy of our 
analytical solution (see Section l4.2.2|) . and also to identify 
terms capable of increasing the accuracy of the solution. 

Table 6: Numerical values of the coefficients Jq{f)p) in the case of 
Janus. The "— " symbol indicates that the corresponding value is 
lower than 1. X 10~®. 



p 




Ji(/3p) 




J3(/3p) 




J5(/3p) 


1 


0, 


.92488 


0, 


.26626 


0.03732 


0.00346 


0.00024 


0.00001 


3 


0, 


.99899 


0, 


.03173 


0.00050 








5 


0, 


.99990 


0, 


.01009 


0.00005 








7 


0, 


.99998 


0, 


.00464 


0.00001 








9 


0, 


.99999 


0, 


.00254 










11 


1 


.00000 


0, 


.00153 










13 


1 


.00000 


0, 


.00099 










15 


1 


.00000 


0, 


.00067 










17 


1 


.00000 


0, 


.00047 










19 


1 


.00000 


0, 


.00034 










21 


1 


.00000 


0, 


.00025 










23 


1 


.00000 


0, 


.00018 










25 


1 


.00000 


0, 


.00014 










27 


1 


.00000 


0, 


.00011 










29 


1 


.00000 


0, 


.00008 










31 


1 


.00000 


0, 


.00007 










33 


1 


.00000 


0, 


.00005 










35 


1 


.00000 


0, 


.00004 










37 


1 


.00000 


0, 


.00003 










39 


1 


.00000 


0, 


.00003 











p 


Jo(fiv) 




MM 


./3(/3p) 






1 





.22708 


0. 


57708 


0.35159 


0.12806 


0.03366 


0.00695 


3 





.98695 


0. 


11367 


0.00652 


0.00025 






5 





.99868 


0. 


.03634 


0.00066 








7 





.99972 


0. 


01673 


0.00014 








9 





.99992 


0. 


00914 


0.00004 








11 





.99997 


0. 


00552 


0.00002 








13 





.99999 


0. 


00356 










15 





.99999 


0. 


00241 










17 


1 


.00000 


0. 


00169 










19 


1 


.00000 


0. 


,00121 










21 


1 


.00000 


0. 


,00089 










23 


1 


.00000 


0. 


,00067 










25 


1 


.00000 


0. 


,00051 










27 


1 


.00000 


0. 


,00039 










29 


1 


.00000 


0. 


,00030 










31 


1 


.00000 


0. 


,00024 










33 


1 


.00000 


0. 


,00019 










35 


1 


.00000 


0. 


,00015 










37 


1 


.00000 


0. 


,00012 










39 


1 


.00000 


0. 


,00010 











B Table of notations 



Table 8: Definition of main notations used in the paper 







Constant of the third Kcplcr law including J2 (sec Formula \b\ ) 


Jl 




Oblatencss of Saturn 


rrii 




Mass of the body where the subscript i characterized the body 


Ms 




Mass of Saturn 


C, 




reduced mass {rrii/ {rrij + m^;)) 


Rz 




Radius of the body where the subscript i characterized the body 


[I] 




the normalized inertia tensor. 


A,B,C 




Normalized (mfi?^) moments of inertia of the whole satellites 
{A< B < C). 


9 




Rotation angle defined as the angle between the orientation 
of the principal axes A and the line of node of the orbit. 


7 




physical libration, oscillation around a mean uniform motion. 


I 




Mean Anomaly 


f 




True Anomaly 


X 




mean longitude of the orbit 


K 




relative mean longitude (Aj — Xe) 


/3, 




Amplitude of the expanding scries of A 


V 




draconic true longitude of the orbit (/ + cj) 


U) 




argument of pcriccntcr 






longitude of pcriccntcr 


n 




longitude of ascending node 


a 




semi-major axis. 


a 




mean semi-major axis 


aj. 




relative semi-major axis (a,/ — as) 


(Xi 




Amplitude of the expanding scries of a 


e 




eccentricity. 


e 




mean eccentricity 


9 




frequency of the longitude of pcriccntcr (tn). 


s 




frequency of the node (SI). 


n 




mean mean motion (A) 


n 




n - g 


u 




orbital libration frequency. 


(7 




frequency of the free libration. 






mean frequency of the free libration 



17 



